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Abstract 

We consider self-similar solutions describing intermittent bursts in shell models of turbulence, 
and study their relationship with blowup phenomena in continuous hydrodynamic models. First, we 
show that these solutions are very close to self-similar solution for the Fourier transformed inviscid 
Burgers equation corresponding to shock formation from smooth initial data. Then, the result is 
generalized to hyperbolic conservation laws in one space dimension describing compressible flows. 
It is shown that the renormalized wave profile tends to a universal function, which is independent 
both of initial conditions and of a specific form of the conservation law. This phenomenon can be 
viewed as a new manifestation of the renormalization group theory. Finally, we discuss possibilities 
for application of the developed theory for detecting and describing a blowup in incompressible 
flows. 

PACS numbers: 47.40, 05.10.Cc 



Estrada Dona Castorina 110, 22460-320 Rio de Janeiro, RJ, Brazil. Phone: +55 212529 5070, 
Fax: +55 212529 5075, E-mail: alexei@impa.br 



I. INTRODUCTION 



The theory of developed turbulence deals with unstable fluid flows at high Reynolds num- 
bers. This theory describes the energy transport from larger to smaller scales until the flow 
is smoothed at small scales due to viscosity. A number of important phenomena are still not 
well understood and described, e.g., the intermittency and anomalous scaling in the inertial 
range l|, |2j . Intermittent bursts, the spatially and temporarily localized fluctuations of the 
velocity field, are principal ingredients of these phenomena. Such intermittent intense events 
may be related to the blowup phenomenon in inviscid flow 0, 4], i.e., to a singularity devel- 
oping in finite time from smooth initial conditions of finite energy in the case of vanishing 
viscosity. We know that, in compressible inviscid flows, the blowup may lead to formation 
of a shock wave. The question of existence of the blowup in incompressible inviscid flows 
(incompressible Euler equations) represents a long standing and still controversial problem, 
see, e.g., [s T | - 

This paper is focused on the study of scaling properties of blowups in hydrodynamic 
flows. Our study is motivated by the success of the shell models of turbulence to describe 
properties of the developed turbulence in the inertial range 8|. The shell models are simple 
"toy" models, where physical space is represented by a discrete sequence of scales taken 
in geometric progression, and the velocity field for each scale is represented by a single 
complex number. Turbulent bursts in shell models are described by self- similar asymptotic 



solutions, which were discovered and described numerically in [9Hl2j. These solutions rep- 
resent the blowup in inviscid shell models, i.e., a singularity develo ping in finite time, which 



we verify using the Beale-Kato-Majda type criterion proposed in 



13] . Note that the self- 



similar solutions are important for explaining anomalous scaling of developed turbulence 
Q, 3- Analogous self-similar solutions were found in more sophisticated cascade models 



of incompressible Euler equations 16 ]. 

It is natural to relate self-similar blowup structures in discrete shell models to blowup 
phenomena in realistic (continuous) hydrodynamic models. In this paper, we accomplish 
this task for a blowup in compressible flows corresponding to the well-known process of 
shock formation. First, we do that for the simplest nonlinear wave equation, the inviscid 
Burgers equation, and then generalize the results for hyperbolic conservation laws in one 
space dimension. A numerical method for detecting self-similar structures in the Fourier 



space is developed by generalizing the approach used in [10j for shell models of turbulence. 
This method works for systems with quadratic nonlinearities, thus, encompassing a large 
variety of hydrodynamic models with quadratic convective terms. 

Self-similar structures in physical space are explained using the renormalization group 
approach. We show that, near the blowup, there exists a limiting wave profile under renor- 
malization of time, space and dependent variable. This limiting function (found in jl^] for 
the case of Burgers equation) turns out to be universal, i.e., independent of initial conditions 
as well as of a particular form of wave equation. Scaling properties of this universal solution 
agree with the Kolmogorov-Obukhov power law, which was not the case in the shell mod- 
els. We compare the obtained results with more sophisticated universality phenomena like 



18 



191 ]; see also 20] for the review of 



critical phenomena of second-order phase transitions 
self-similar singularities in problems of fluid dynamics. 

The paper is organized as follows. In Section 2, we describe self-similar blowup solu- 
tions in shell models of turbulence. Section 3 describes generalization of these results for 
the Fourier transformed Burgers equation. In Sections 4 and 5 we return to the physical 
space and explain universality of renormalized blowup solutions for compressible flows. Sec- 
tion 6 presents the renormalization group approach. Section 7 discusses application of the 
developed theory for incompressible flows. Section 8 summarizes the contribution, and the 
Appendix contains some technical derivations. 



II. SELF-SIMILAR STRUCTURES IN SHELL MODELS OF TURBULENCE 

Shell models of turbulence are simple dynamical systems that demonstrate important 
properties of developed turbulence in the Navier-Stokes equations at high Reynolds numbers. 
In these models, the turbulent velocity field with wave numbers k in the spherical shell 
k n < ||k|| < k n+1 {n = 1,2, . . .) is represented by a single complex number u n called shell 
velocity. The shell wave numbers are chosen as k n = ko\ n with the most popular value 
A = 2 for shell spacing. Equation for the shell velocity is 

du n /dt = Q n [u] - fk 2 n U n + /„, (1) 

where Q n [u] describes the quadratic nonlinearity, v is the viscosity, and /„ is the external 
forcing applied to the first shells. It is assumed that only several neighboring shells interact, 
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as described by the nonlinear term of the form 

Q n [u] = i{k n+1 u n+2 u* n+1 - ek n u n + x u* n _ x + (1 - £)/c n _ 1 w n _iW n _ 2 ) 



(2) 



corresponding to the Sabra model [21]. For the GOY model 22. 23 1. the nonlinear term is 
obtained by taking all shell speeds in fl2]) with an asterisk and substituting (1 — e) by (e — 1). 
The parameter is often chosen to be e = 0.5, in which case both the "energy" ^ \u n \ 2 and 
the "helicity" £(-l) n fc„| conserved in the system with no viscosity and no forcing, 

v = f n = 0, reflecting basic properties of the Navier-Stokes equations. 

Fig. [TJ shows typical behavior of the shell speeds. It is characterized by a sequence of 
intermittent bursts, which have a self-similar structure shown in the inset. In each burst, 
speed fluctuations propagate to smaller scales (shells with larger n), until the viscosity 
becomes dominant and smooths the oscillations. Self-similar structure of such bursts was 
in |a- 



described 



11| by considering equations with vanishing viscosity v = 0. It was shown 



(3) 



numerically that the shell velocities behave asymptotically for large n as 

u n (t) = -u c (t c - tf'^F (u c (t c - t) 1/z k n ) , t<t c . 

We presented this expression in the different (equivalent) form, which is convenient for 
comparison with continuous models below. In ([3]), t c is the final time and u c > is the 
scaling parameter. The function F(k) is purely imaginary and defined up to scaling F(k) h-> 
u c F(u c k), which reflects the symmetry of the shell model (JTJ, (J2J) with v = f n = under the 
scaling k n h- > k n /u c , u n \- > u c u n . This function is universal for all the bursts. Expression (J3]) 
describes identical structures, which appear in subsequent shells renormalized by scaling of 
time and speeds. Numerical computations provide the value 



0.719 



(4) 



for the GOY model with e = 0.5 



10| . For the Sabra model, the best fit of rescaled velocity 



profiles for a specific turbulent burst gave z = 0.75±0.02 [11] . However, as we will see below, 
the exponent z has the same value (jlj) for both the GOY and Sabra models. This exponent 
is different from z = 2/3 corresponding to the Kolmogorov-Obukhov law u n ~ k n derived 
on dimensional grounds. This difference has implications for explaining anomalous scaling 
in developed turbulence [ll |. 

At the final time t c , the value of shell velocity u n is finite for any n and, thus, no singularity 
is created for a particular shell. Despite of that we can show that expression (J3|) describes a 




FIG. 1. Temporal dynamics of the Sabra shell model with 20 shells, e = 0.5, ko = 2 -4 , v = 2 x 10~ 6 , 
fi = 0.05. Curves with smaller values of \u n \ (in average) correspond to larger n. The inset shows 
the self-similar structure of the shell velocities with n = 7, ... ,13 for the largest turbulent burst. 

blowup, i.e., a singularity developing at t c for the infinite series of shells 1 < n < oo. This 
can be done using the Beale-Kato-Majda type criterion of blowup in inviscid shell models, 
which requires [13j 

lim / sup k n \u n {s)\ds = oo. (5) 

t— >tc Jo n 

Using (j3J), we express 

K\u n \ = {t c -t)- l y\F{y)\, y = u c (t c -t) 1/z k n . (6) 

For small t c — t and k n = koX n , one can always choose (large) n such that 1 < y < A. Hence, 

sup k n \u n (t) | > a(t c - ty 1 , a = min y\F(y)\, (7) 

where a is a positive constant. This expression implies divergence of the integral in ([5]). 

It is interesting to find a type of solution, which appears at t = t c . Using ([3]) and 
expressing \u n \ similarly to ([6]), we find 



lim \u n \ = u z c k z n lim y z F(y). 

Since each u n is finite at t c , the limit in ([8]) is finite and \u n (t c 
constant b. It follows that, at the blowup time t c , the norm 

t 1/2 

\U\d = I J2 k2 n\ U ^ 2 
\n=l 



(8) 

bk^ 1 for some positive 



(9) 



is infinite for d > 1 — z = 0.281 and finite for d < 1 — z. Thus, we checked directly the 



blowup condition, which was introduced in 



13] as the divergence of the norm with d > 1. 



Moreover, since 1 — z < 1/3, the (week) solution at t > t c does not satisfy the energy 
conservation criterion ljj. Hence, the blowup time t c may be viewed as a starting moment 
of the energy- cascade. 

Expression (J2J) can be explained using dynamical scaling of time and shell speeds. This 
theory was proposed in [10!], and it is reviewed in the rest of this section. In our presentation, 
we changed the form of derivations and final self-similar expressions facilitating development 
of the analogous theory for continuous hydrodynamic models. 

Let us write equation (JTJ) with v = f n = for the scaled speeds as 

dv n /dt = N n [v), v n = ik n u n , (10) 

with the nonlinearity 

N n [v) = -\~ 2 v n+2 v* n+l + evn+iv^ + (1 - e)\ 2 v n -iv n -2- (11) 



For phenomena far from the first shells, i.e., for large n, we can consider equation (fTTj) 
defined on the infinite lattice n 6 Z. Then, system (fTUl) . ffTTj) is invariant under translations 
n i — y n+1. For simplicity, we will assume that all v n are real, so that u n are purely imaginary. 
In this case, the GOY and Sabra models are equivalent. 

A translation-invariant system possessing a positive quadratic invariant can have traveling 
wave solutions. Such a system can be obtained from (TTOj) by introducing the dynamically 
scaled speeds w n (r) with the new time r as 

t = t c - J exp ^- jf A(r")dr'' \ dr', v n = exp Qf A(r / )rfr'^ io n , (12) 

where the blowup time t c corresponds to r — > oo. It is straightforward to check that 

dw n /dr = N n [w] - Aw n , (13) 

where iV n [w] has the form ffTTj) written in terms of w n instead of v n . Equation ([TBI conserves 
the sum w 2 if we take 

A(T) = 5> n Ar>]/X>n- ( 14 ) 

n n 

Note that transformation f|T2|) - ffT^|) is valid for any quadratic nonlinearity N n . 



FIG. 2. Solid lines correspond to the shell model (Sabra or GOY) and dashed lines correspond to 
the Burgers equation. Shown are the profiles W(£) of solitary waves (|15p and (|22p . imaginary part 
of the universal functions F(k) in ([3]) and (j23j) . and the asymptotic form ([3]) of shell speeds u n (t) 
near the blowup with t c and u c adjusted to match the inset in Fig. [TJ 



Numerical computations 10[ showed that equations f ll~3|) . ffTTl) have an asymptotically 
stable solution in the form of a solitary wave moving toward large n (small scales) with 
constant speed s, i.e., 

w n (r) = W((n- sr) log A) (15) 

for all t and n; the function VF(£) vanishes as £ — > ±oo. The extra factor log A is introduced 
here, which will be convenient for comparison with continuous models below. Solution ([15]) 
written in term of the original shell speeds u n (t) yields the self-similar expression ([3]) with 
the function 

F(k) = ik^exp ( [ A(r')dr'] W(-sr log A), z = — ?— [ ' A{r)dr, (16) 
\Jo J log A J 

where r is related to t = t c — k z implicitly by (TT2|) ; see Appendix for the derivation. The 
function W(£) of the solitary wave ffTBT) and the corresponding purely imaginary function 
F(k) given by f|T6|) are presented in Fig |2j The graph on the right of Fig |2] shows shell 
velocities ()3]) describing the asymptotic form of the blowup presented in the inset of Fig. [TJ 



III. SELF-SIMILARITY OF BLOWUP FOR BURGERS EQUATION 

Self-similar structure of the blowup should have an analogue in continuous hydro- 
dynamic models. First, it is natural to analyze simple models having well-known blowup 
structures. In this section, we start with the Burgers equation u t + uu x = vu xx + /, where 
v is the small viscosity, / is the external forcing, and subscripts denote partial derivatives. 



The Burgers equation models various phenomena in compressible fluids, multiphase flows in 
porous media, etc. For the Fourier transformed wave profile u(k) = {2tt)~ 1 J u(x)e~ tkx dx, 
this equation takes the form 



du 
dt 



/oo 
(k — p)u(p)u(k — p)dp — k 2 vu + /. (17) 
-oo 



For real functions u{x) decaying sufficiently fast at infinity, the Fourier transform u(k) is 
regular and u(—k) = u*(k). Thus, we can consider only positive k from now on. Additionally, 
we assume that u(k) is purely imaginary, which means that the function u(—x) = —u(x) is 
odd in physical space (we return to the general case in the next section). 

In order to study a blowup, we take v — f — 0. Then, there is a global invariant J u 2 dx, 
which corresponds to j \u\ 2 dk in the Fourier space. (Also, one can treat J uu x dx = as a 
second global invariant analogous to helicity.) Hence, equation (|T7|) can be related to the 
shell model ([I]), fl2]) with the same reasoning as for the Navier-Stokes equations. The shell 
wave numbers k n = ko\ n growing in geometric progression induce the expression k = in 
the continuous formulation. For the logarithmic coordinate £, we write (I17p with v — f = 
as 

Ov/dt = N[v], v(0 = -ik 2 u(k), (18) 
with the real function v(£) and the nonlinear term 



N[v] = W 

<T=± J - 



30 £2 

v{v)v{()— d V, 
PQ 



(19) 



(20) 



k = e^, p — ae v , q = = \k — p\, 

where we used the relations 

u{p)dp = u(p)\p\drj = u(\p\)pdr) = (i/p)v(rj)dr}, 

(k - p)u(k - p) = qu(q) = (if q)v(Q 

following from u(—k) = u*(k) = —u(k) for the purely imaginary u(k). One can see the 
analogy of equations f lT8|) . ( IT9"]) with the discrete system f JTOj) . (ITT]). In particular, N[v] is 
quadratic with respect to v and invariant under translations + A^. 

We can write (fTBl . ( Tl9l) in the form similar to ( Tl3l) . ( [141) as 

/oo / /-oo \ —1 

wN[w]d£ ( / w 2 di j . (21) 
-oo \J — oo / 



The real function w(t,£) is related to v(t,£) by expressions ffT2l) . where the subscript n must 
be dropped. One can check that equation f l2~T|) conserves the integral J w 2 d^. Hence, it 
allows traveling wave solutions. 

Numerical integration of ( 1211) shows that there exists a stable solution 

w(t,0 = W(£-st), (22) 

representing a solitary wave moving with constant speed s in positive direction of £-axis, 
i.e., toward small scales x ~ 1/k = e~^. The profile W(£) of this wave is shown in Fig. [2] 
by the dashed line, and one can notice that it is very close to the corresponding profile ( IT5|) 
obtained for the shell model of turbulence and shown by the solid line. Derivations given in 
the Appendix lead to the self-similar solution analogous to (j3j) of the form 

u(t,k) = u 2 c (t c -tf z - l F(u c (t c -t) l ' z k), t<t c , (23) 

with arbitrary u c > and 

F(k) = ^w(logk + ^y z=± (24) 

where the constant A is given by (121j) . The exact Kolmogorov value of the exponent 

z = 2/3 (25) 

is found in this case. A graph of the purely imaginary function F(k) is presented in Fig. [2] 
(dashed line). 

We see that the traveling wave solution ffT5l) for the shell model of turbulence has a direct 
analog (|22|) for the inviscid Burgers equation, with the corresponding self-similar solutions 
(EJ), (1TB]) and ( 1231) . ( 1241) . The result ( 1251) suggests that the anomalous value of z in (jlj) may 
be attributed to the approximate nature of the shell models. Physical interpretation of the 
obtained self-similar solutions is given in the next section. 

IV. SELF-SIMILAR BLOWUP STRUCTURE IN PHYSICAL SPACE 

The self-similar blowup solution just described must correspond to the blowup (formation 
of a shock wave) in the inviscid Burgers equation 

u t + uu x = 0. (26) 

9 



The solution u(t, x) constructed by the method of characteristics has the implicit form 



x = x + u (x )(t - t ), u = u (x ), (27) 

where u(t ,x) = Uq(x) is a smooth initial condition and Xq is an auxiliary variable. For 
spatial derivative of u(t, x), we have 

du _ du/dxp _ u' (x ) 

dx dx/dx 1 + u' (x )(t — t ) ' 

The denominator vanishes at t = to — 1/u (xq). This yields the well-known result that 
the classical solution blows up along the characteristic with the minimum negative value of 
u' (x ), followed by the formation of a shock wave. 

Let us choose the origin of time and space so that the blowup singularity appears at 
t — x — 0, and consider a smooth solution in the interval to < t < 0. Also, we can take 
u = at the singularity, which can be achieved by the transformation Uo(0)t and 

it i — y u -\- Mo(0), which leaves (|2T)|) invariant. In this case, the initial condition Uo(x) satisfies 
the conditions 

uo(0) = 0, *o = l/«d(0) < 0, «o'(0) = 0, <(0)>0, (29) 

which ensure that u' (x) has a negative minimum at the origin. Note that, for the odd 
functions u(x) considered in the previous section, the conditions uq(0) = u'q(0) = are 
satisfied automatically. 

For small Xq, we use (|29|) to expand 

uo{xo)to = xo + ^J^x 3 + o{x 3 ). (30) 



Using this expansion in ([27]) yields 

6u'(0) 



u"'(0) 

x = ut Xq + o(xq), u = u' (0)xo + o(x ). (31) 



Equivalent ly, 



x = ut-cu 3 + o(m 3 ), c = ^tt^L > 0- (32) 



6(m'(0)) 4 

Consider now the renormalization operator Q\ acting as 



u x (t, x) = G x u(t, x) = A 1/3 w(A- 2/3 t, X^x). (33) 
10 



It is easy to see that u\ is a new (renormalized) solution of (|26|) . Thus, (|33|) represents 
a symmetry of (126]) . Multiplying both sides of (132]) by A and making the substitution 
1 1 — y A~ 2 / 3 t, x i — V \~ 1 x yields the equation for u\(t,x) as 

x = u x t-cul + \o((\- 1/3 u x ) 3 ) ■ (34) 



The last (correction) term contains powers A 1 n ^ 3 with n > 3. Hence, it vanishes in the 

nn 

limit of large A, and (I34"j) takes the exact form [17|, |2Cj 



X = M^t - CM 3 ^ (35) 

for the limiting function 

Uooit^x) = lim u\(t,x). (36) 

A— >oo 

Clearly, the function Uoo(t, x) is self-similar, 

£a«oo = Moo. (37) 

We write solution of fl35|) as 

Uoo(t,a;) = u c (-t) 1 / 2 F(u; 1 (-t)- 3 / 2 a;), u c = c" 1 / 2 , £ < 0, (38) 

where F(x) is defined uniquely by the equation 

x = —F — F 3 , (39) 

as one can check by substituting (|38|) into (|35|) . Note that the factor w c reflects the scaling 
symmetry of the Burgers equation (126]) . 

We see that, near the blowup, the renormalized solution u\(t,x) is described asymptoti- 
cally by the function F(x) in (13"81 . Fourier transform of (13"B1 yields the self-similar expression 
(T23|) with z = 2/3 and t c = 0, where u(t,k) is the Fourier transformed function Uoo(t,x). 
This explains self-similarity of blowup in Fourier transformed Burgers equation with the 
exact Kolmogorov value z = 2/3 of the scaling exponent. Also, the function F(k) describing 
self-similarity in the Fourier space is the Fourier transform of F(x) defined in (1391) . This 
completes analytical description of the self-similar blowup structure for the inviscid Burgers 
equation. 

For connection with the classical theory, we can consider the limit t — > with fixed 
x 7^ 0, when the argument of the function F in fl38|) becomes large. Equation fl39l) yields 

11 
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FIG. 3. (a) Convergence of the renormalized profiles u\(to,x) (black curves) to the universal form 
Uoo(to,x) (bold red curve) given by (f38]h (f39|) for the inviscid Burgers equation ([26]) . The initial 
condition is uq{x) = — ( v / 7r/2)erf x at to = — 1> and A = 1, 10, 10 2 , 10 3 . (b) Similar convergence 
result for blowup of a simple wave in ideal polytropic gas. The bold red curve represents the 
universal limiting function Loo (to, %) given by (|44|) . (|38|) . (|39p and back curves are the renormalized 
profiles U\(to,x) = Q\U(to,x) corresponding to the gas density variation U. 



F(x) ~ —x 1 / 3 for the large argument x 3> 1- Using this asymptotic formula in (|38l) . 
we obtain the cubic root relation ~ — ul^x 1 ^ 3 at the singular point, which is well- 
known. Note that (!38l) represents a much stronger result, since it suggests that the whole 
renormalized profile (with derivatives of any order) becomes universal. 

Fig. [3^l presents a numerical example of convergence of the renormalized solution u\(to, x) 
to the self-similar form fl38l) with the increase of A for t = — 1 and the initial condition 
Uq(x) = — (y/n/2)erf x. We also confirmed numerically that F(k) is the Fourier transform 
of F(x). 

We conclude that the self-similar solution ( 1231) obtained by the Dombre-Gilson method 
successively capture the universality of renormalized profiles near the blowup in the inviscid 
Burgers equation. This universality result turns out to be much more general, as we will see 
below. 
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V. UNIVERSALITY OF BLOWUP IN COMPRESSIBLE FLOWS 



Compressible (or multi-phase) inviscid flows in one space dimension are modeled by 
hyperbolic systems of conservation laws. In scalar case, U Gl, we consider the equation 

dU_ df(U) _ 
dt dx 

with a smooth flux function f(U) and initial condition Uq(x) at t = t . A classical solution 
of equation fj40|) has the form 

x = x + f(U (x ))(t-t ), U = U (x ), (41) 

which reduces to solution ( 1271) of the Burgers equation by the substitution u = f'(U). 
We assume that the coordinates for time and space are chosen such that (1291) holds for 
uq(x) = f'(Uo(x)), which implies the blowup at t = x = u = 0. Assuming that u = f'(U) 
is locally invertible as U = g(u) with g(Q) = and g'(0) = l//"(0), we have the blowup at 
{7 = for the original state variable. We will use the relation 

g x u n (t, x) = \ 1/3 u n (\- 2/ X X-'x) = A (1 -" )/3 <(t, x) (42) 

with the renormalization operator Q\ given in (1331) . Expanding g(u) in Taylor series 

1 u 1 

= p(0) + ^(0)« + -g"(0)u 2 + ... = — + -g"(0)u 2 + ■ • • , (43) 



Ueofax) = lim Q x U(t,x) = lim ^(«(t,ar)) = ^(t, x)/f(0), (44) 

A— >oo A— >-oo 



we find 



where all terms Q\u n with n > 1 vanish for A — > cxd due to (I42p . 

Formulas (1331) . (|44|) with the function u^it.x) given by (|38|) . (1391) demonstrate strong 
universal character of the blowup in scalar ID conservation laws. We see that, when the 
singular point is approached, the renormalized wave profile U(t, x) takes the universal form 
independent both of the initial condition Uq(x) and of the flux function f(U). In particular, 
derivatives of all orders turn out to be universal for the renormalized solution near the 
singular point. Note that expressions (l32l)- ([36l) relate the blowup with a cusp catastrophe 



in the space (t, x, u), see 24j. Therefore, universality of the blowup just described represents 



a property of this catastrophe. 
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As an example, let us consider formation of a shock in a simple wave solution for one- 
dimensional flow of ideal polytropic gas. The gas density p(t, x) in this wave is described 
implicitly by 

X-X 7 + 1 r-r- ( 7 -l)/2/ s / s , A -s 

7 ^7Po (x Q ), p = po{x ), (45) 



t - 1 7-1 



where po(x) is the initial condition at t = t , see, e.g., [25]. We will use the values 7 = 5/3, 



A = 3/5, and po(x) = 2 — arctanx. Then expressions (|4"5]) take the form 



x = xq + 4(2 — arctanxo) 1 ^ 3 (t — t 



(46) 



p = 2 — arctanxo, 

which can be written in the form (1411) for the density variation in moving frame 

U (t, x) = pit, x + xx + vxit- t )) - pi (47) 

with 



/'([/) =4(tf + p 1 ) 1 /3_ Vlj 
C/ (x) = 2 — arctan(x + xi) — p\. 
The quantities t , xi, p\ and v x are chosen such that f'(0) = and the function 



(48) 



u {x) = f'{U {x)) = 4(2 - arctan(x + Xi)) 1/3 - v 1 (49) 

satisfies conditions (|29|) corresponding to blowup at t = x = with [7 = 0. The values of 
these quantities are 

t = -1.155, x x = 0.183, pi = 1.818, ^ = 4.882. (50) 

Convergence (l4*4"l) of the renormalized profile of gas density variation U\(t, x) = Q\U (t, x) to 
the universal limiting form is demonstrated in Fig. [3b. 



VI. RENORMALIZATION GROUP 



The scaling universality of blowup can also be explained using the renormalization group 
approach. For this purpose, we consider Q\ in ( 1331 as an operator acting in the space 
of solutions of the Burgers equation ( |26l) with initial conditions satisfying ( 1291) . Then the 
system evolution can be seen as the action of Q\ combined with renormalization of space x 
and state u. The operator Q\ defines a differentiable one-parameter group with the property 
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<?Ai<?a 2 — G\!+x 2 - The self-similar solution Uoo(t, x) in (!38|) . ( )39|) represents a stationary point 
f )37|) of the renormalization group operator. This stationary point (more precisely, a set of 
stationary points u^it^) determined up to a scaling constant u c ) is asymptotically stable 
in the sense of Lyapunov for A — > oo considered as "time" This stability condition 

implies the universal limit (156]) . 



In the case of a general conservation law the renormalized function U\(t,x) = 

G\U(t,x) is a solution for a conservation law with the renormalized flux function f\(U) = 
A 2 / 3 /(A~ 1 / 3 f/). Thus, Q\ defines a one-parameter group acting in the functional space of 
solutions and fluxes as (U, f) (U\,f\). The universality f H4|) of the blowup is explained 
by the fact that Q\ has the asymptotically stable stationary point (U, f) = (Moo(t,x), U 2 /2) 
corresponding to solution ( |38l) of Burgers equation ( 1261) . 



The role of the operator Q\ is similar to renormalization group operators in other physical 
theories. For the inner scale, x ~ ut ~ (— 1) 3//2 with small £, the dynamics is governed by the 
universal function, which is a stationary point of Q\. This is analogous, e.g., to the stationary 
point of the renormalization group operator, which determines critical exponents in second- 
order phase transitions |l8] ■ At larger spatial scales, there is no universality and the solution 
depends on the initial condition Uq(x) as well as on the flux function f(u). This is analogous, 
in turn, to the phenomenological Landau theory of second-order phase transitions valid at 
larger (though still small) deviations of temperature from a critical value 26 1. 



Up to now, we considered compressible flows, where the blowup corresponds to the well- 
known phenomenon of shock formation. Could the presented approach be used for the study 
of incompressible flows? We address this question in the next section. Such an application 
is strongly motivated by the conclusion, which was made in jj] based on the scaling analysis 
of a blowup in numerical solutions of incompressible Euler equations: "it could be that the 
antiparallel configuration is universal as t — > T and r < R. For large r there is a helical 
regime that might not be universal..." This hypothesis is exactly what the renormalization 
group theory described above suggests for the blowup solution. 
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VII. DOMBRE GILSON SCHEME FOR INCOMPRESSIBLE EULER EQUA- 
TIONS 



The procedure described in Section 3 works for evolutionary systems with quadratic 
nonlinearities written in the Fourier space. This covers a large class of hydrodynamic models 
with quadratic convection terms. Let us consider the incompressible Euler equations written 
for the Fourier transformed velocity u(£, k) as 

du/dt = -iP J d 3 p (u(p) ■ k) u(k - p) (51) 

with the incompressibility condition k • u = 0. Here P = I — kk T jk 2 with k = ||k|| is the 
projector onto transverse vector fields k • u = 0, and the incompressibility condition is used 
to remove the pressure and u(p) • p terms. Considering logarithmic spherical coordinates 
£ = log A;, Ofc = k/k, we write equation for the scaled speed in the form similar to (ITS)) as 

dv/dt = N[v], v(£,o fc ) = -^ 4 u(k), (52) 
N[v] = P J d 2 o p d V ^ (v(p) • o„) v(k - p), (53) 

where p = ||p||, q = ||k — p||, r\ — logp, o p = p/p and the integration is carried over 
— oo < f] < oo and the unit sphere o p G S 2 . One can check that this equation is invariant 
under translations £ H- £ + A£ of the vector field v(£, which represent scaling in the 
Fourier space, k i-> e A ^k. We assume that v is a real vector field corresponding to a flow with 
the symmetry u(— x) = — u(x). Since the nonlinearity N[v] is quadratic, we can transform 
()53p to the form similar to (1211) as 



dw/dr = N[w] - Aw, A =(J d 2 o k di ||w|| 2 ^ J d 2 o k d^w ■ N[w] 



(54) 



where r and w are given by (jl2|) with v n and w n substituted by v and w. System (jo! 
is translational-invariant along direction £ and preserves the norm J d^o^d^ ||w|| 2 . Hence, 
it admits traveling wave solutions of the form w(r, £, Ok) = W(£ — sr, Ofe). Such solutions 
determine a self-similar blowup described in physical space by an expression analogous to 



which was considered first by Leray j27j] in this context. 
Numerical computations for system (1541) were performed, where the vector spherical har- 
monics expansion with coefficients depending on £ was used for the function w(£, o^). Such 
expansion has advantages for numerical computation of the operator iV[w]. Details of the 
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numerical method and its realization will be published elsewhere. Here we comment only on 
the conclusions related to the presented approach. First of all, no traveling wave solutions 
were found. On the other hand, we observed that at least two characteristic scales describe 
formation of a singularity, since the solution was concentrated in a small region on the unit 
sphere o^. £ S 2 with increasing time r and wavenumber k. This conclusion agrees with the 
numerical analysis in (4] suggesting the two scales, r ~ (t c — t) and (t c — t) 1 ^ 2 , where t c is 
the time of anticipated collapse. 

Note that equations (lo"^|) account for self-similar solutions characterized by a single scale 



x ~ k~ x = e~^. Such solutions were ruled out in 28] assuming sufficiently fast decay of 
the vorticity at infinity in physical space. We see that the Dombre-Gilson scheme must be 
generalized in order to capture several scales. Possible generalizations are discussed in a 
separate paper [29 ]. 



VIII. SUMMARY 

Representation the blowup in the shell models by a traveling wave solution in the log k 
space suggested by Dombre and Gilson [10| is shown to be physically relevant for continuous 
hydrodynamic models. For compressible inviscid flows, a similar wave describes the universal 
scaling properties of shock formation. 

We showed that solutions of hyperbolic conservation laws describing compressible inviscid 
flows in one space dimension have universal scaling structure when approaching a finite-time 
singularity. The limiting renormalized wave profile is described by a universal function in 
physical space, or by a solitary wave moving with constant speed in logarithmic coordinates 
of the Fourier space. This universal function is independent both of initial conditions and of 
specific expressions for fluxes of conserved quantities. We also explained this phenomenon 
using the renormalization group theory. Finally, we showed that the Dombre-Gilson scheme 
can be used for the study of blowup in incompressible inviscid flows. 
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APPENDIX 



Using ([TOD, ([I2D and ([ISD, we fin d 



u n {t) = -ik~ l v n {t) = -ik~ l exp (j A^drj W((n - sr) log A), (55) 

Note that expression (I14p with w n from f lT5|) is invariant with respect to the shift n t- > n + 1, 
t i— > t + 1/s. Hence, the function A(r) is periodic with period 1/s, and we can write ( 155]) as 

u n (t) = -k~ 1 \ nz F 1 (r-n/s). (56) 

Here the exponent z with the function Fi(r) are defined as 

1 f 1/s 



log A 



A(r)dT, (57) 



o 



F 2 (t-t c ) = F 1 (T)=iexp(J A(t') dr^j W{-sr log A), (58) 

where the function F 2 will be used later, and t is related to r by (|12j) . 

Let us find the time t = t n corresponding to r — n/s. Using (jl2j) . we obtain 

exp f - y A(r")dr" j dr ' = ~ J exp ( ~ / A(r")rfr"j dr' 

= -A n2 exp f- A(r")rfr"J rfr = A n2 (t - t c ), 

(59) 

where we used the change of variables r' h- )■ r' + n/s and periodicity of A(r"). According to 
(EHD and (J59J), we have F x (r - n/s) = F 2 (t n - t c ) = F 2 (X nz (t -t c )). Using this function in 
( 1561) with /c n = fc A n and taking u c = k$ , we derive the expression 



( (t) = -K z k z - X F 2 {ko Z K(t - t c )) = -u c {t c - tf' z - l K z - x F 2 {-K z ) , (60) 



where K = u c (t c — t) l l z k n . This expression with F 2 from (|58j) and 2; from (1571) reduce to 
(UJ) with r related to t = U - K z by (PT2D. 
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In the continuous case, we substitute k n by k = (so that A = e). We have A = const 
in (|2ip for the traveling wave solution (122]) . In this case we obtain expressions analogous to 
dSSD (jnOD in the form 

u(t,k) = ik~ 2 v(t,0 = k z - 2 F 1 (r-Z/s) = k z ~ 2 F 2 {(t - t c )k z ) (61) 

with 

z = A/s, F 2 (t - t c ) = Fi(r) = ie AT W(-sr). (62) 
For nonzero constant A, the first expression in 0121) yields the following relations 

e~ Ar log(A(t e - 1)) 
f = t c -_ , r = ^. (63) 

We reduce expression (1BT]) to the form (B31 with m c = 1 by introducing the function 

F(k) = k z ~ 2 F 2 {-k z ). (64) 

Note that an arbitrary factor u c > in (|23|) reflects the scaling symmetry of the Fourier 
transformed Burgers equation (fTT|) with u — f — 0. Using fl63l) in fl62l) yields 

Expression ([24]) for F(k) follows from (]64]). (1651). 
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